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ON ■ We review the Skyrme model approach which treats baryons as solitons of an ef- 

fective meson theory. We start out with a historical introduction and a concise 
Ph! discussion of the original two flavor Skyrme model and its interpretation. Then 

we develop the theme, motivated by the large Nc approximation of QCD, that the 
^ ■ effective Lagrangian of QCD is in fact one which contains just mesons of all spins. 

When this Lagrangian is (at least approximately) determined from the meson sector 
it should then yield a zero parameter description of the baryons. We next discuss 
the concept of chiral symmetry and the technology involved in handling the three 
flavor extension of the model at the collective level. This material is used to discuss 
properties of the light baryons based on three flavor meson Lagrangians containing 
just pseudoscalars and also pseudoscalars plus vectors. The improvements obtained 
by including vectors are exemplified in the treatment of the proton spin puzzle. 
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1. Historical background and motivation 

The Skyrme model was born around 1960 in a series of increasingly more detailed papers Q. 
At that time the prevailing dynamical model of nuclear forces was that of Yukawa which had 
been formulated in the 1930 's. Still to come was the concept of fractionally charged quarks and 
much further in the future was the recognition that the correct theory of strong interactions 
binds these quarks together with non-Abelian (color) gauge fields. 

In the Yukawa theory, of course, the nucleons are introduced as fundamental fermion fields 
while the spin zero pion fields are postulated to provide the "glue" which binds protons and 
neutrons into nuclei. This model is acknowledged to work reasonably well as a description of 
the long range interactions of nucleons and the prediction of the existence of pions has been 
amply confirmed. 

Skyrme's innovation was to provide a model in which the fundamental fields consisted of just 
the pions. The nucleon was then obtained, in the initial approximation, as a certain classical 
configuration of the pion fields. The seeming contradiction of making fermi fields out of bose 
fields was avoided by arranging the classical field configuration to possess a non-zero "winding 
number". In modern language this "Skyrmion" is an example of a topological soliton. Such 
objects are solutions to the classical field equations with localized energy density ||. They play 
an important role nowadays in many areas of physics and the papers of Skyrme are justifiably 
recognized as pioneering milestones in this development. 

The years following this original idea saw the particle physics community actively investi- 
gating the approaches of quark models, flavor symmetry, current algebra, chiral dynamics, dual 
resonance models and finally color gauge theory to the problem of strong dynamics. Evidently 
Skyrme's model was lost in the rush. However the novelty of the model did stimulate a few 
interesting papers || f|, ||, || before the more recent wave of activity in the area. 

At first glance, it might appear that a Lagrangian model built out of only pion fields could 
not be more different as a description of the nucleons from the current picture of three "valence" 
quarks containing a trivalent "color" index and bound together through their interaction with 
SU(3) gauge fields. Remarkably, it has turned out that the Skyrme model is in fact a plausible 
approximation to this QCD picture. This may be understood as follows. 

In QCD the gauge coupling constant has an effective strength which decreases for interac- 
tions at high energy scales (asymptotic freedom) but which increases at the low energy scales 
which are relevant when one considers the binding of quarks into nucleons and other hadrons. 
Thus the application of standard perturbation theory techniques to the problem of low energy 
interactions is not expected to be reliable and in fact has not produced definitive results. A 
natural alternative approach which retains the possibility of using perturbation theory is to 
imagine that the strong underlying gauge couplings bind the quarks into particles which may 
possibly interact with each other relatively weakly. At low energies these particles should evi- 
dently comprise the pseudoscalar meson fields (pions when restricted to two "flavors"). Then 
it is necessary to formulate some effective Lagrangian for the pions. Certainly the Lagrangian 
should be restricted by the correct symmetries of the underlying gauge theory. These must 
include an (approximate) SU(Nf) flavor symmetry p|], where Nf is the number of light flavors. 

But there is another symmetry which plays a crucial role. At about the same time that 
Skyrme was contemplating the model under discussion the correct formulation || of the struc- 
ture of the effective weak (beta-decay etc.) interaction was discovered. It was noted that this 
interaction treated the left and right handed components of fermions on a completely separate 
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basis. If this distinction is maintained at the level of the strong interactions one should impose 
a "chiral" left handed SU(3) flavor x right handed SU(3) flavor symmetry on the effective low 
energy Lagrangian of mesons. A consequence of this larger symmetry is that the meson multi- 
plets must contain scalar as well as the desired low-lying pseudoscalar particles. This seemed a 
bit of an embarrassment in that the pseudoscalars are very light while the scalars were not well 
established and presumably heavy. Hence the possibility of a degenerate symmetry multiplet 
is implausible. Nevertheless it was realized || that the situation was likely to be similar to 
that met in the BCS theory of superconductivity in which the vacuum (ground state) is en- 
ergetically favored to exist in a non-symmetric state. This "spontaneous breakdown" picture 
predicts, in the absence of a needed small explicit symmetry breaker, exactly zero mass for the 
pseudoscalars at the same time that the scalars are massive. In fact it may be formulated using 
a so called "non-linear realization of chiral symmetry" in such a way that the scalars do not 
appear at all The prototype Lagrangian density for this picture is 



f 2 

c=iftT(d IM u& i u*) + ..., (1.1; 
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where U = exp(\/2i0// 7r ), <\> being the 3x3 matrix of the ordinary pseudoscalar mesons and 
f n = 93MeV the "pion decay constant". U is a unitary matrix which transforms "linearly" 
under the chiral transformations. Possible higher derivative and symmetry breaking terms 
have not been explicitly written here. It was demonstrated a long time ago that just this term 
compactly summarizes the low energy scattering of pseudoscalar mesons JO]] . Improvements to 



this term form the basis of the "chiral perturbation theory approach" [p^] . Now it is believed 
that a picture like this is expected from fundamental QCD. However the same Lagrangian was 
earlier written by Skyrme (in the two rather than three flavor case) in order to explain the 
nucleons Q before the present justifications for it were known. 

Even in the framework of the chiral Lagrangian given above it would seem that there is no 
special a priori reason not to explicitly add baryons in a chiral symmetric manner rather than to 
build them out of the mesons. Indeed there have been many papers over very many years which 



do just this with reasonable phenomenological results [[y|. Nevertheless there is an indication 
from fundamental QCD that the soliton treatment of the baryon is more natural. This arises 
from an attempt [[L4| to consider 1/Nc, the inverse of the number of colors in the gauge theory, 
as a possible expansion parameter for QCD which might be meaningful even at low energies. 
In this approach the product g' 2 = g 2 N c , where g is the gauge theory coupling constant, is held 
constant, 't Hooft |14| showed that for large Nc QCD may be considered as a theory of mesons 



weakly interacting in the sense that scattering amplitudes or quadrilinear coupling constants 
are of order g e g = 1/Nc- Since the baryon mass must start out proportional to Nc (noting that 
the baryon in the Nc model is made of Nc quarks) it means that the predicted expressions for 
baryon masses should start out as the inverse of this coupling constant. In the framework of 
the (non-relativistic) mean field treatment Witten jn| not only pointed out that the baryon 
masses indeed grow linearly with Nc but also that baryon radii and meson-baryon scattering 
amplitudes are of the order N c while baryon-baryon scattering is of O(Nc)- He in particular 
recognized that this inverse behavior with g e g is just the usual signal that the baryon state in 
question is a soliton of the effective meson theory. 

Naturally, one wonders how these "modern" justifications for the Skyrme approach relate 
to Skyrme's original motivations. We are fortunate in having available a reconstructed talk on 
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just this topic by Skyrme fl6 |. He mentioned three motivations: 1) The idea of unifying bosons 



and fermions in a common framework. 2) The feeling that point particles are inconsistent 
in the sense that their quantum field theory formulation introduces infinities which are only 
"swept under the rug" by the renormalization process. 3) The desire to eliminate fermions from 
a fundamental formulation since fermions have no simple classical analog. What seems more 
fascinating is his awareness that there were probably some "hidden" influences pushing him 
toward the soliton picture. Directly, these came from his fascination with Kelvin's idea that the 
various atoms should correspond to vortices of different connectivities in some underlying liquid. 
In turn, his interest in Kelvin was sparked at an early age by the presence of a tide prediction 
machine, designed by Kelvin and constructed by his great-grandfather, still occupying space 
in his great-grandfather's house. An interesting account of this aspect is given in a paper of 
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Thus it seems that Skyrme's motivations were not those currently used to justify his model. 
In particular it appears that he did not choose his Lagrangian model to describe spontaneous 
breakdown of chiral symmetry. Rather the non-linear form was chosen to insure that the pions 
were "angular" variables which would give multi-valued functions; the crossing of different 
sheets of these functions might then correspond to singularities which would realize the baryons. 
The evident "moral" of this historical discussion is just that interesting ideas have an uncanny 
way of turning out to be useful and true. In this spirit, we would like to continue with the 
application of Skyrme's ideas to current research on baryon physics, making use of current 
motivations but trying to avoid getting enslaved by them. 

2. The Skyrme model for two flavors 

In this section we will present the basic technology of the Skyrme model for baryons. The 
starting point for the construction of a soliton solution is the non-linear sigma model Lagrangian 
( |1.1| ) already introduced in the previous section. As we require a finite energy density the chiral 
field U must approach a constant value at spatial infinity. We are free to choose this to be 
unity, i.e., 

U(r,t)-+1. (2.1) 

This can be considered a mapping from compactified coordinate space, a three-sphere S 3 , to 
the space which is described by the unitary, unimodular matrix U, namely SU(Nf), where 
Nf denotes the number of flavors. In the case of two flavors the target space is isomorphic to 
S 3 . The mappings S 3 — ► S 3 fall into distinct equivalence classes. This signals the existence of 
soliton configurations because members of different classes cannot be continuously transformed 
into one-another. The equivalence classes are characterized by the winding number. This 
number counts the coverings of the target space and is the charge / d 3 xB associated with the 
topological current 



B » 247r 2 efMUpa tr 



(U ] 3 U U) (U ] 8 P U) (U ] d a u)] . (2.2 



When later discussing the three flavor case we will see that this topological current indeed 
equals the baryon number current. 

Although these topological considerations allow the existence of soliton solutions it turns 
out that the dynamics of ( |1.1| ) do not lead to static stable classical solutions. This can be 
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deduced from a simple consideration, known as Derrick's theorem [18]. Assume Uo(r) to be 
such a solution. The static energy of U (Xr), obtained from the Hamiltonian of (1.1), would 
then be 



E^ ) [U (Xr)} = J Er > [U (r) 



-,(nl CT )r 



(2.3) 



which does not have a minimum at A = 1, in contradiction to the assumption. In order to 
obtain stable solitons Skyrme added a term to the Lagrangian which is of fourth order in the 
derivatives, 



£(Sk) 



1 



32e 2 



tr 



UtdpU, Wd u u] [u ] d^U, U ] d u u] ) 



(2.4) 



Here e is the dimensionless "Skyrme constant" . Although this term is quartic in the derivatives 
it was arranged to be at most quadratic in the time-derivatives. This makes the quantization 
feasible. It is now apparent that a scaled configuration may well lead to a minimum of the 
energy functional 



E^mXr)} = jE^'lUoir)] + XE^ K> [U (r)] 



,(111(7) r 



,(Sk)r 



(2.5) 



at A = 1 provided the configuration Uq(v) satisfies E^ l<7> [Uo(r)] = E^ k \Uo(r)]. 

A priori the Euler-Lagrange equations of motion for the chiral field Uo{r) are highly non- 
linear partial differential equations. To simplify these equations Skyrme adopted the famous 
hedgehog ansatz 



Uo(r) = exp (zr • rF(r)) 



(2.6) 



where r represents the Pauli matrices. This form may actually be traced back to the old 
"strong coupling" theory |TJ| . Upon substitution of this ansatz the energy functional turns into 
a simple integral involving only the radial function F(r), 



E[F] 



dx 



Henceforth this radial function will be called the 
chiral angle. In eq ( |2.7| ) a prime indicates a deriva- 
tive with respect to the dimensionless coordinate 
x = ef n r. In this manner we have completely 
extracted the dependence on the model parame- 
ters. Imposing F(oo)=0 and noting that / d 3 rE>o = 
(F(0) — F(oo))/ti leads to the boundary condition 
F(0) = 7r for a unit baryon number configuration. 
The profile function depicted in Fig. 2.1 minimizes 
(|2.7|) and is obtained numerically. The energy ob- 



(x 2 F 12 + 2sin 2 F) + sin 2 F (2F' 2 + ^f^j j 



(2.7) 




Figure 2.1: Chiral angle 



tained by substituting this solution into (|2.7|) is 
found to be E = 23.2nf n /e @. 

As the ansatz (|2.6| ) is not invariant under separate spatial or flavor rotations this field 
configuration does not yet describe states with good spin and flavor quantum numbers. As 
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the first step towards generating such states, time-dependent collective coordinates A(t) are 
introduced which describe the spin and flavor orientation of the hedgehog, 

U(r, t) = A(t)U (r )A*(t), A(t) e SU(N f ) . (2.8) 

Note that the hedgehog structure causes rotations in coordinate and flavor space to be equiv- 
alent. For generality we assumed an arbitrary number of flavors. In the two flavor case this 
configuration yields the Lagrange function 

L(A,A) = ^a 2 n 2 -E cl (2.9) 

where the quantity ft measures the time dependence of the collective coordinates 

A\t)jA{t)= l -r-n. (2.10) 

The constant of proportionality, a 2 is computed as a spatial integral over the chiral angle to 
be a 2 = 53.3/(e 3 / 7r ). Computing the spin as the Noether charge of spatial rotations yields 
J = a 2 Q. Apparently the system displays all the features of a rigid top. In that language O 
and a 2 are denoted as the angular velocity and the moment of inertia, respectively. 

In the second step the collective coordinates are elevated to quantum variables. Again this 
is completely analogous to the quantization of the rigid top and gives the quantization rule 
[Ji, Jj] = ieijkJk- For the hedgehog ansatz in SU(2) spin and isospin are related via the adjoint 
representation of the collective coordinates, i.e. Ii = —DijJj with Dij = (l/2)tr(riArjA^) due 
to the equivalence of the respective rotations. Hence only states which have identical spin and 
isospin are allowed in the spectrum. These are the nucleon (J = J = 1/2) and the A-resonance 
(/ = J = 3/2). Finally the Hamiltonian for the collective coordinates is given by 



H coU = E d + -^J 2 = E d + —I 2 (2.11} 



1 _o _ 1 
'coll — •E'cl + 7. , -. 

2cr 2a 
which yields the A-nucleon mass difference 



M A -M N = — . (2.12) 

Using the physical value for the pion decay constant, f n = 93MeV requires us to choose 
e ~ 4.75 to reproduce the empirical mass difference of 293MeV. Substituting e ~ 4.75 into eq 
( |2.7| ) yields the classical nucleon energy E = 23.2nf n /e ~ 1430MeV. This is not in especially 
good agreement with the experimental value of about 939MeV. However the following points 
must be kept in mind: 

(i) The meson Lagrangian consisting of eq Q1.1Q plus eq ( |2.4| ) contains only pseudoscalars. We 
would expect that other low mass mesons (notably the vectors) should also be included. 
The large Nc expansion |]T4 , |i~5H requires an infinite number but common sense suggests 



a reasonable approximation for explaining hadronic physics up to about lGeV would 
be to keep those mesons with masses up to this value. Certainly the predictions in the 
mesonic sector of the theory are noticeably improved by the inclusion of vector mesons. 
The consistency of the overall picture requires accurate predictions both in the mesonic 
and baryonic sectors of the effective theory. 
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(ii) In nature there are three rather than two "light" flavors and this aspect should be included 
in a realistic formulation. (This feature also makes more transparent the origin of the 
topological current eq (|2.2|) .) Furthermore the effects of flavor and chiral symmetry 
breaking mediated by the finite values of the quark masses have not yet been taken 
into account. 

(iii) Order of N c corrections to the nucleon mass which have the structure of the Casimir 
effect in field theory have also not been included. These quantum contributions to the 
energy have been estimated to be negative and of the order of a few hundred MeV, 
predicting a total nucleon mass at the order of the experimental value [E0|. Nevertheless 



one should be cautious about these quantum corrections, after all the Skyrme model is 
not renormalizable, leaving a logarithmic scale dependence of the "renormalized" Casimir 
energy. It seems that at best the quantum corrections can be computed in a scenario 
compatible with the chiral expansion. 

We will postpone the discussion of a variety of nucleon (and other baryons') properties until 
after we have treated the more general case of flavor 577(3). 

To end this section on the basics of the Skyrme model we would like to briefly discuss the 
consistency of the Skyrme model with the large Nc picture of QCD. In section 1 we have 
already noted that the quadrilinear coupling between mesons scales like 1/Nc- To check this 
behavior it is convenient to expand the non-linear a model Lagrangian ( |1.1| ) in powers of the 
pion field: 

i<VT • 3"tT + {(7T ■ d.Tzf - TT^TT • C^Tt} + O (V) . (2.13) 

Since the quadrilinear coupling constant is l// 2 we deduce that f T ~ y/Nc- This agrees with 
general arguments fl5fl . Similarly the Skyrme term (2A) provides a quartic pion interaction 



with the coupling constant l/(e 2 /^) which implies e ~ 1/y/Nc- Hence the classical energy ( |2.7| ) 
grows linearly with the number of colors as asserted from the corresponding generalization of 
QCD. Moreover, without flavor symmetry breaking large Nc QCD predicts the baryons of 
different J to be degenerate [51]. This is perfectly consistent with the mass formula ( |2.11| ) 



because the moment of inertia also grows linearly with Nc as is indicated after eq ( |2.10|) , hence 
the second term in ( [2.1 1|) behaves like 1/Nc- 



Actually, the understanding of the Nc expansion for baryons involves some subtleties. Con- 
sider the construction of large Nc baryons in the quark model. The lowest lying baryons are 
made of Nc (taken to be odd) quarks in a totally antisymmetric (i.e. singlet) color spin state 
with no orbital angular momentum. One expects particles of all total angular momenta from 
J — 1/2 to J — Nc/2 to be obtained. In agreement with the spectrum of eq (|2.11|) we expect 



I = J for these particles and an infinite number of them as Nc — > oo. The trouble is that there 
is no experimental evidence for any I = J = 5/2, 7/2 etc. particles. This may be interpreted 
as evidence that Nc = 3 in nature. Still, the large Nc expansion is useful if one computes a 
quantity which exists in the Nc = 3 theory as a (presumably quickly convergent) Taylor series 
in 1/Nc- For the specific case of the higher excitations J = 5/2, 7/2, . . . the above treatment 
of the rotational modes seems inadequate because the rotational energy gets as large as the 
classical contribution. By including these modes in the Euler-Lagrange equations the widths 
of these higher excitations have been estimated to be comparable to their masses p2j. This 
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makes a particle interpretation of these states problematic suggesting that they are artifacts of 
the collective quantization method employed rather than of physical relevance. Possible caveats 



for these calculations are the instability of these configurations against emitting pions [23] and 
that the results are only obtained by analytical continuation in the spin variable. 

3. Chiral symmetry and its breaking 

In this section we will briefly discuss the concept of chiral symmetry which represents a 
guiding principle for extending the Skyrme model. Attention will be limited to those aspects 
of this large subject which have direct relevance to the study of Skyrmions. The basic idea is 
to construct a model of meson fields which "mocks up" as many symmetries and properties of 
the fundamental QCD Lagrangian as possible. 

3.1 The QCD Lagrangian 

Let us first recall the matter piece of the QCD Lagrangian 

N f 

/=i 

N f 

= J2 {Qf,L + 94) Qf,L + ?/,R + 94) Qf,R - m f (9/,L?/,R + Qf,RQf,h)} ■ (3- 1 ) 

/=1 

Here is the matrix representation of the gluon fields and g the quark-gluon coupling. Most 
notably we have introduced the chiral representation for the QCD current quarks 

?l,r=^(1T7 5 )? (3-2) 

of each flavor. 

Strictly speaking, the quark mass terms are not part of the QCD Lagrangian but arise 
from the Yukawa terms of the full microscopic theory of nature. A major unsolved problem 
is to understand the resulting pattern of quark mass parameters. The phenomenologically 



determined masses [fPJ] are m u ~ 5MeV, ~ 9MeV, m s « 120 — 170MeV, m c ~ 1.5GeV, m& ~ 
4.5GeV and m t « 175GeV. This random looking perturbation of the "strong interaction" 
plays a crucial role in determining the nature of elementary particle physics. In the region up 
to about lGeV it is not possible to produce particles containing c, b or t quarks. Then it is 
usually a good approximation to simply drop them from the theory. Other approximations 
are useful when dealing with the subspace carrying the flavor quantum number of a single 
"heavy" quark p5| . Furthermore in the sector of the three "light" quarks u, d, s it turns out to 



be fundamental to neglect the u, d, s masses as a first approximation and include their effects 



as a perturbation [^6|. This is reasonable because the light masses are less than the quantity 
Aqcd ~ 250MeV, the scale below which the QCD effective coupling gets extremely large. 

In the case m/ = the Lagrangian (|3.1|) specialized to the three light quarks has the global 
chiral symmetry 

q u 

t4(3) x C/ R (3) : q L — > Lq L and q R — > Rq R with Ql,r = \ Qd \ ( 3 -3) 



Qa 



L,R 
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where L and R are each 3x3 unitary matrices. Using Noether's theorem on the classical 
Lagrangian then yields the conservation of the eighteen vector and axial vector currents 



Jij = Qjl^Qi and j§ >B = ^7 M 75%, (3.4) 

where the latin indices run over u, d, s. These currents play an important role in the theory of 
weak interactions. 

Now a major discovery of quantum field theory is that consequences of the classical field 
equations of motion (which can be used to verify the conservation of the Noether currents) do 
not necessarily hold at the quantum level. It is necessary to consider whether there exists a 
suitable regularization of the divergent diagrams of the theory which maintains the classical 
relations. In the present case, the axial singlet current[|, j\ = Ql^i^ /2)<Z is not conserved 
even for massless quarks. Rather its divergence is proportional to the gluon field tensor times 
its dual. This is a result of the well-known Adler-Bell-Jackiw (ABJ) triangle anomaly P? | 



contained in the loop-diagrams shown in figure 3.1 




Figure 3.1: Adler-Bell-Jackiw anomaly. The arrows indicate quark lines, the curly lines refer 
to the gauge bosons and the dashed lines denote the coupling of the axial singlet current J5 . 

The net result is that the true global symmetry of the massless quantum theory is not 
C/l(3) x Z7r(3) but Uy(l) x S77 L (3) x SUr(3). The singlet vector symmetry, Uy(l) corresponds 
to baryon number conservation. 

A similar situation emerges when external c-number flavor gauge fields are added to the 
massless QCD Lagrangian in order to further probe its structure. Then the so-called non- 
Abelian anomaly yields non-zero covariant derivatives of the S77l(3) x S77r,(3) currents propor- 
tional to certain combinations of the corresponding external gauge fields . The non-Abelian 
anomaly will be noted to have important consequences for the theory of Skyrmions. 

Although the true symmetry of massless three flavor QCD is S77l(3) x S77r.(3) x Uy(l), the 
resulting symmetry of the physical states of the theory is further reduced to SUy(3) x Uy(l) 
by the "spontaneous breakdown" mechanism. In this mechanism the vacuum state is not 
invariant under the full symmetry group. The massless QCD vacuum is characterized by a non- 
vanishing "condensate" (q u q u + qdQd + q s q s ) 7^ 0. Under an infinitesimal chiral transformation 

1 Here A a , a = 1, . . . , N^ — l denote the Gell-Mann matrices of SU(Nf) while A refers to the singlet generator. 
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TV 2 — 1 N 2 — 1 

L = 1 + iJ2 a =o £L^ a /2 and R = 1 + iZ) a =o e R^ a /2 the variation of this quark-bilinear is 
found to be 

S (qq) =H e l- 4) ?Ly?R - gRy<?L )=(<£- 4) <?y*7 5 g • (3.5) 



Clearly the condensate is invariant only for the subgroup L = R in eq ( |3.3D which is a vector type 
transformation. This explains the physical SUy(3) x Uy(l) invariance. Note that the right hand 



side of eq (|3.5| ) represents pseudoscalar objects. These are "zero mode" fluctuations of the above 
vacuum configuration and the corresponding massless pseudoscalar particles are designated 
Nambu-Goldstone bosons. Their scalar chiral "partners" - which would be degenerate in mass 
were it not for the spontaneous symmetry breakdown - are not constrained to be massless. 
This splitting of low-lying pseudoscalars and scalars expected from massless QCD seems in 
qualitative agreement with the experimental situation. 

The SUy(3) x Uy(l) invariance (so-called "eightfold way") of massless QCD is, of course, 
further broken when the effects of non-zero quark mass terms are included. For later purposes 
it is convenient to rewrite the quark mass terms as: 



C^ = -^±^qMq with M= 2 -^l + yX 3 + 1 -^X 8 . (3.6) 



Here characteristic quark mass ratios are defined by 

2m s m u - m d 
x = , y = . (3.7) 

m u + m d m u + m d 

In the limit y = 0, the theory possesses SUy(2) or isospin invariance. 

3.2 Effective Lagrangian of pseudoscalars 

The simplest way to mock up low energy QCD is to employ a 3 x 3 matrix field Mij which 
transforms under the chiral group in the same way as the bilinear quark combination qjRqu- 
It has the decomposition M = S + iP into hermitian scalar and pseudoscalar components. A 
chirally invariant Lagrangian is 

C = -tr^M^M 1 ") - V(M, M r ), (3.8) 

where the potential V is a function of invariants like tr(MM*), \v( M M' M M : ) etc. Spontaneous 
breakdown to SUy(3) is implemented by choosing V to have a minimum such that (M) = 
const, x X Then the scalar fields S become massive and can be "integrated out" by imposing a 



chirally invariant constraint |]10[ . This represents a transition from the linear to the non-linear 
sigma model. Formally we may use the "polar decomposition" of the matrix M into unitary 
and hermitian factors M = HU . Setting H — > const, x 1 then results in eq ( |1 . 1| ) again. 

In principle, a scenario of this sort can be derived by adding a term like q^Mq^ + h.c. to 
the QCD Lagrangian and then integrating out the quark fields^. As a result one is left with a 

2 In the case of QCD this seems to be impractical. However, the simpler Nambu-Jona-Lasinio jj) model 
for the quark flavor dynamics nicely exemplifies how a meson functional can be constructed by integrating out 
quark degrees of freedom p3| . 
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complicated action functional for M and, by further eliminating H as above, for U. Note that 
U inherits the chiral transformation property of the quark bilinear: 

U — ► LUB) . (3.9) 

It is this identification of the transformation properties which provides the important link of 
the effective chiral theory to QCD since it in particular implies that the (Noether) currents 
must be identified. In turn, matrix elements of these currents will yield the static properties of 
hadrons. 

Of course, in the chiral limit we demand the effective meson theory to be strictly invariant 
under the transformation (|3.9|) . Since WU = 1 only derivative terms can appear and the 
leading one is just the non-linear a model ( |1 . 1| ) . Clearly also the Skyrme term fl2.4|) is invariant 
under (p|). 

At this point one essential ingredient is still missing to ensure that the chiral field, U = 
exp(i$) describes pseudoscalar fields, $. This requirement demands the parity transformation 

parity , 

U-^rf. (3.10) 

However, it is straightforward to verify that the Skyrme model Lagrangian is invariant under 
( |3.10| ) and r — > —r separately. On the level of the equations of motion we can easily break 



this unwanted extra symmetry by adding a term which contains the Levi-Civita tensor. With 



a 



{d^U)U^ we write 



yfya" + ... + 5Xe flupa a p a u a p a (7 = , (3.11) 



where the ellipsis refers to the contributions from the Skyrme term ( J2.4Q which have the same 
symmetries as those from the non-linear a term ( |1 . 1| ) . Unfortunately, the additional term 
cannot be easily incorporated in the effective Lagrangian since it does not correspond to the 



variation of a local term. Witten suggested []30| to include it at the level of the action because 
the variation of 

T = \f tYe lxvpaT a tl a v a p a <T a T (3.12) 

and the use of Stoke's theorem yields the desired term in the equation of motion (|3.11|) provided 
the boundary of the five dimensional manifold, M 5 is taken to be Minkowski space, i.e. dM 5 = 
M4. The choice of M5 is not unique because its complement has the same boundary. In order 
to nevertheless have a unique action the constant A = 2 ~J" 2 , n e Z must be quantized^. It 
is interesting to study the physical relevance of ( |3.12| ). Expanding in the meson fields $ and 
employing again Stoke's theorem reveals that it describes processes with at least five different 
pseudoscalarsQ like K + K~ — > 7r + 7r°7r~. As such processes were first discussed by Wess and 
Zumino J31] who essentially found a power series expression for (ft. 12 ), the term is commonly 
named after themQ. 

3 The reader is referred to the literature j3(| to see the analogy to Dirac's quantization of the magnetic 
monopole. 



4 For that reason the term (3.12) vanishes in the case of two flavors which only has four different pseudoscalar 
fields. 

5 Some authors refer to the Wess-Zumino term in its gauged form i.e. with external sources. Unless otherwise 
noted we will always understand the Wess-Zumino term to be ( |3.12 ). 
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There are further important consequences of the Wess-Zumino term ( |3.12|) which can be 
read off after generalizing it so its variation with respect to external (electro-weak) gauge 
transformations [3C, 32|] yields the non-Abelian anomaly [p8|| . After appropriately including 



the corresponding gauge boson fields two striking features are observed: 

(i) A contact interaction for the decay 7r° — > 77 is contained in the gauged Wess-Zumino 
action. On the quark level this process is described by the ABJ anomaly involving the 
diagrams of figure |3.1| with the external lines representing photons. Identifying that 
result with the Wess-Zumino term requires setting n = Nq, i.e. the Wess-Zumino term 
is proportional to the number of colors. 

(ii) The linear coupling to the Uy(l) gauge boson represents the baryon number current. 
Indeed it turns out that this current is identical to the topological current in eq 



The mocking up of the effects of the Ux(l) anomaly in QCD involves the SU(3) singlet 
pseudoscalar particle rj and will be discussed later when we treat the proton spin puzzle. 

We must also take account of the effects of the finite quark mass terms ( |3.6| ). These trans- 
form according to the chiral SU^(3) x SUr(3) representation: 3 x 3* + 3* x 3. Note that 
the matrix M. = M.^ (neglecting the possibility of strong CP violation) may be considered a 
"spurion" for this transformation property. Then the minimal symmetry breaking piece of the 
effective Lagrangian reads 

£sb = tr [M [-(3' (d^Ud^U + U ] d^Ud^U ] ) +5' (U + U ] - 2)] } , (3.13) 

where f3' and 8' are two numerical parameters. The 5' term is required to split the pseudoscalar 
meson masses while the f3' term is required to split the pseudoscalar "decay constants". The 
decay constants f a are defined from the axial vector matrix elements {0\j^\(pa,p) = ifaPfi 
Working in the isospin invariant limit, the parameters /?', 5' and x can then be extracted 



from the knowledge of meson properties [33 



4 „, 4 ( f K Y . 4 



ml = -5', m l = —5\l + x) and ^ = 1 + —(3'{1 - x) . (3.14) 

Jtt JK V JTT / Jtt 

This represents the essential input when discussing the Skyrme model for three flavors. To sum 
up, the Lagrangian of only pseudoscalars which we shall use for discussing Skyrmions consists 
of the sum of Q, 0, (gig ) and (|3~T3D . 



In the chiral perturbation theory approach |12| essentially the most general chirally invariant 
Lagrangian is written down and ordered in powers of dd ~ M.. For example, the leading terms 
are eq ( |1.1| ) and the 5' term of eq ( |3.13| ). The next-to-leading terms include: 

[trid^Ud^)} 2 , ti(d fl Ud l/ U j )tT(d^Ud u U j ) } 

tv(d lx Ud^d 1/ Ud v U^), ti(d^Ud^U ] )ti(M{U + f/ f )), ti(d^Ud^U ] M{U + U r )), 
[ti{M{U + U^)] 2 , [tv(M(U - U ] )f MMU ] MU ] + MUMU), (3.15) 

each with its own coupling constant. Note that the combinations of terms which can not be 
manipulated (by use of various matrix identities) to become a single trace are suppressed in the 
l/Nc expansion. This procedure also entails absorbing the divergent parts of loop corrections 
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in the coefficients of the listed terms. The result is a joint power series in energy and the quark 
masses which can be expected to be very accurate quite near the tttt threshold in the case of 
pion-pion scattering for example. But going higher in energy is very difficult in this scheme. 
Furthermore it seems that many of the coefficients mainly simulate the low energy effects of 
vector meson exchanges. We shall also work with a meson Lagrangian which includes the vector 
particles directly. This may be thought of as a start on the approach of constructing the leading 
(Born) term of the 1/Nc expansion, which should include mesons of all spins. 

3.3 Effective Lagrangian of pseudoscalars and vectors 

We will follow the so-called massive Yang-Mills approach jS^] for introducing the vector 
meson nonet into the Lagrangian of pseudoscalars in a chirally invariant manner. In this 
approach both vector and axial vector fields are formally introduced as gauge fields (yielding 
invariance under local chiral transformations). Then globally invariant mass-type terms which 
break the local chiral invariance are included. Finally, as in the transition from the linear to 
the non-linear sigma model discussed in section 3.1 above, the (heavier) axial vector mesons 
are eliminated by a chirally invariant constraint. 

We introduce two multiplets with spin one, and A^ which we demand to transform 
under ( jOD as left- and right-handed fields, respectively, 



— > L UJ + -aj Lt and A* 



r K + 



(3.16) 



This allows us to define a covariant derivative for the chiral field and field tensors, 



F L,R 



dM-igA h u U + igUA* 



d,A u 



L,R 



d u A 



L,R 



,iL,R /iL,R 



(3.17) 
(3.18) 



which transform homogeneously under 
number of derivatives read 



3|). The chirally invariant terms with a minimal 



tr 



(D^UyD^U 



tr 



and tr 



F L JJF K,IJ,V U^ 



In addition we can have mass-type terms for the vector mesons 



tr 



A\A^ + A|M R " 



and tr 



A h UA K ^U ] ' 



(3.19) 



(3.20) 



which are still invariant under global chiral transformations. Of course, many more terms with 
higher derivatives could be written down at the expense of more undetermined parameters. 
Now, it is our aim to construct an effective model for the vector mesons only; at present we 
are not interested in the axial-vector mesons. We have to find a mechanism to eliminate the 
latter without violating the chiral symmetry. This can be accomplished by choosing a special 
"gauge" for the vector fields A^' R , 




(3.21) 



Here p M is a matrix field with Nj components. For example, in the case of two flavors it includes 
both the p and u> mesons via p M = p • t + u^X In the case of three flavors, this matrix field 



13 



is supplemented by the K* and cj) mesons. Most importantly we have introduced the "square 
root", £ of the chiral field, U = which yields the chirally invariant relation 




A L = U{A* + ~dJtf. (3.22) 



It is actually this so-called unitary constraint which eliminates the axial-vector fields in favor 
of the vector fields p without spoiling chiral symmetry. 

It is interesting to study the behavior of the p meson under chiral transformations. To start 
off, we recognize that the transformation of £ introduces the matrix K which is defined by [[34 1 



i — ► = K^B) . (3.23) 



Clearly this leaves the transformation law of the chiral field (|3.9| ) unchanged. Note that in 
general the matrix K is a position dependent quantity because of £. Demanding now the 
symmetry transformation 

(3.24) 

causes the fields y4 L,R to transform exactly like left- and right-handed vector fields. 

Within the unitary gauge the various terms listed above in (|3.19 ) and (|3.20|) are no longer 




independent. Introducing the homogeneously transforming combinations 

P» = d£$ + tL%t and ^ = ^ + ^-(0^-^) (3.25) 

the terms up to two derivatives can be combined to a chirally invariant Lagrangian of vectors 
(and pseudoscalars) 



£vm — tr 



(3.26) 



where we have used the fact that the coefficient of the term quadratic in the p meson field 
is related to the vector meson mass m p = 770MeV. Upon expanding the square-root field 
£ in powers of the pseudoscalar field, one finds that the Lagrangian ( |3.26|) contains the prni 
coupling, 



m 2 p 



£p™=2^2P M -(7rx3"7r) . ( 3 - 27 ) 

which can be utilized to fix the coupling constant g ~ 5.6 from the known decay-width of the 
process p — > rnv. 

Terms which involve the Levi-Civita tensor e^ upa are also of great interest for the Skyrme 
model. For their presentation it is most useful to introduce the notation of differential forms: 
A R = A^dx^, d = d p dx^, etc. . Since the left- and right-handed "gauge fields" are related 
via the unitary constraint Q3.22 ) the number of linearly independent terms, which transform 
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properly under chiral transformation as well as parity and charge conjugation, is quite limited 



A L a 3 , dA L aA L - A L adA L + A L aA L a , 2 (A l ) 3 a + -A L aA L a . 

v ' 9 



(3.28) 



For convenience we have again made use of a M = (dnU)W = £p£L Of course, including these 
terms in the model Lagrangian will introduce three more parameters: 71,72 and 73. A suitable 
presentation of this part of the action is given in terms of p and R (employing again the notation 
of differential forms) 



wz 



+ 



tr 

m 4 V6 



7i + 



Rp*- g l2 F(p)\pR-Rp}-g 2 [ l2 + 2 l3 }R i p) , (3.29) 



where T-wz is given in fl3.12|) . In ref [35] two of the three unknown constants, 71,2,3 were 
determined from purely strong interaction processes like lu — > 3tt. Defining h = -2^2^/3, 
9vv<t> = 912 and k = 73/72 the central values h = ±0.4 and g~vv<f> — ±1-9 were found. Within 
experimental uncertainties (stemming from the errors in the uj — <fi mixing angle) these may vary 
in the range h = —0.15, . . . , 0.7 and gvv<j> — 1-3, . . . , 2.2 subject to the condition \gw<j>~ h\ ~ 1-5- 
The third parameter, k could not be fixed in the meson sector. From studies |3(| of nucleon 
properties in the two flavor model it was argued that k ~ 1 represents a reasonable choice. 

The sum of the space integral of (|3.26|) and (|3.29|) comprise the chirally invariant part of 
the effective action we shall use for discussing the soliton in the vector meson model. Note that 
the second piece of (|3.29j ) can be gauged with external fields |35|] so as to make no contribution 
to the non-Abelian anomaly. The first piece T^z then correctly supplies the non-Abelian 



anomaly. Furthermore the second piece of ( 3.29 ) stabilizes the soliton without the need for 
including the Skyrme term Q2.4p . 

We must still include the effects of symmetry breaking due to finite quark masses in the 
vector meson system. To leading order in the symmetry breaking, an appropriate term which 
behaves properly under chiral transformations, can be constructed by analogy to the last ex- 
pression in (|3.20|) 



- a'tr \M (A^UA r " + A*lPA L i* 



(3.30) 

This leading contribution not only distinguishes between the p and K* masses but also con- 
tributes to the different decay constants of the pseudoscalar mesons via the unitary gauge 
( 3.21 ). The reader may consult ref [33] for recent discussion of higher order symmetry breaking 



terms. 

We would like to end this section on including vector mesons by noting that the same La- 
grangian is obtained within the so-called hidden gauge approach |J7), once the same symmetries 
are required. This shows that these two approaches are in fact identical. 

3.4 Other aspects 

We expect that baryons should appear as solitons of the large Nc effective meson Lagrangian 
for any number of flavors Nf. In the case of three (or more) light flavors the Wess-Zumino 
term guarantees, as discussed in section 3.1, that the baryon number (|2.2j ) is obtained in a 
self-contained manner from the Lagrangian. This can be used to check that the soliton indeed 
has the correct baryon number. 
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Now in the two flavor case, the same kind of soliton solution exists. However the Wess- 
Zumino term vanishes identically so we cannot similarly check its baryon number in a self- 
contained way. The situation is even more peculiar for Nf = 1. There the Skyrme model 
represents a mapping S 3 — > S 1 which does not contain topologically stable configurations. 
However, we are not forced to use an effective Lagrangian of the same form. In this case it is 
probably more realistic to construct the Lagrangian by including isoscalars like the spin-0 o- 
field and the spin-1 cu-field. Such a Lagrangian might have a soliton solution (not necessarily 
topological) but a check of its baryon number may also not be available in a self-contained way. 
These examples seem to indicate that the form of the relevant effective Lagrangian may have 
a non-trivial Nf dependence (at least for small Nf). 

Another interesting question, related in the sense of understanding whether physical features 
of the solitons can be traced to particular pieces of the effective Lagrangian, concerns the 
stabilization of the soliton. In section 2 we noted that the Skyrme term (|2.4j ) was introduced 
precisely for this purpose. There is an often mentioned "derivation" of this term from the 
piece of the vector meson Lagrangian (|3.26|) above which goes as follows. In a large mass 



expansion, m p — > oo the equation of motion for the vector meson field simply becomes R p = 0. 
Substituting this into the remainder of the vector meson Lagrangian (|3.26|) 



yields exactly the Skyrme term ( |2.4j ) with the identification g = e. Although the numbers 5.6 
and 4.75 are in reasonable agreement there is one caveat to this appealing derivation of the 
Skyrme term. While the Skyrme model does yield stable solitons, however, for arbitrary large 
but finite m p the model ( |3.26|) does not contain stable soliton solutions. Thus one seems to 
have achieved stabilization merely by approximating a model in which stabilization does not 
exist. Clearly we have not obtained a "physical origin" for the stabilization mechanism. As 
mentioned in the previous section, the second piece of (|3.29|) stabilizes the soliton in the vector 



meson Lagrangian without a need for the Skyrme term. It is also possible that, as in the case 
of the s-wave ground state of hydrogen, stability is achieved at the quantum, rather than at 
the classical, level. Several investigations of this possibility have been made based on just 



the non-linear sigma model term (|1.1|) , although an assumption on the allowed chiral profiles 
seems to be required. 

4. The Skyrme model with three flavors 

It is well established J7| that the neutron (n) and proton (p) belong to a multiplet with 
six other members (the iso-singlet A, the iso-doublet 5 and the iso-triplet £). To try to 
understand n and p alone is to look at only a small piece of a large picture. Thus we must 
consider the three flavor generalization of the treatment in section 2. First (in the present 
section) we shall consider the Lagrangian of pseudoscalars alone, discussed in section 3.2. The 
new features arise from the inclusion of flavor SU(3) symmetry breaking terms (see ( |3.13| ) 



together with (|3.6| ) and (|3.7|) ) as well as the Wess-Zumino term (|3.12| ). Both of these features 



involve non-trivial extensions of the formalism and interesting "physics" . 

The first step towards including the strangeness degrees of freedom is to actually take the 
chiral field to be a U(l) ® 577(3) matrix. To be precise, the three flavor chiral field is defined 
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as 



U (x) = exp 



V3A 



r/o exp (i$) . 



(4.1) 



While the singlet field 7/0 is separated the matrix field $ now not only contains the pion degrees 
of freedom but also the kaons and the non-singlet component of the 77 fields, 



8 [K 
a=l /t 



V 



7T 



i? 



K + 
K° 



76*18 / 



(4.2) 



Here A a denote the Gell-Mann matrices. Note that in the presence of derivative-type symmetry 
breakers (e.g. the j3' term in (|3.13|) ) the normalization of the fields gets shifted; the "physical" 
fields are gotten by multiplying the fields above by some constants as Z 7r 7i + , Z K K + , etc.; 
similarly the physical decay constants are Z n f 7T = 93MeV, Z^fx ~ 113MeV, etc. For the Z's 
we have 



Jit / 



J K 



4(1+^ 

Jtt J 



etc. 



(4.3) 



We clearly need a suitable generalization of the Skyrme ansatz (|2.6f) . It turns out that it 
is correct to just embed the SU (2) hedgehog in the SU (3) matrix. Flavor symmetry breaking 
implies that field configurations which have non-zero strangeness possess a classical energy 
which (at least in the unit baryon number sector) is larger than that of a zero strangeness 
configuration. Thus we choose the embedding: 



Un(r) 



exp (it ■ rF(r)) 
















(4.4) 



Hence the classical energy will not be modified and the soliton profile, F(r) is that in figure 2.1. 
The effects of the strange degrees of freedom are hence visible when states with baryon quantum 
numbers are generated via the collective coordinate approach. 

The collective coordinate matrix A(t) defined in eq (|2.8|) is now taken from SU(3) and in 
analogy to eq (|2.10| ), now leads to eight angular velocities, 



A ] it)jA(t) = l -j:x a n a 



(4.5) 



In addition to the angular velocities fl a the adjoint representation 

D ab = -ti (\ a AX b A^ 



(4.6) 



of the collective rotations, A(t) will be important, in particular in the context of flavor symmetry 
breaking. 



17 



Substituting U = A(t)Uo(r)A^ (t) into the pseudoscalar Lagrangian of section 3.2 without 
the symmetry breaker gives rise after a spatial integration to the collective Lagrangian 



L Skyrme (A A) + L WZ (A A) = \a 2 + b 2 E^l-^§^s- E cl . (4.7) 

1 i=\ 1 «=4 2 V3 

The SU(2) moment of inertia a 2 remains unchanged while the moment of inertia j3 2 for rotations 
into the strange directions is a new functional of the pseudoscalar fields. The fact that the eighth 
component of the angular velocity vector does not appear quadratically in (PO) is a consequence 
of [Uq, Xs] = 0. The term proportional to BQ$, where B is the baryon number arises from Twz- 



In order to obtain it we make use of the separation [39| 

iN, 



Twz[U] = T WZ [U ] - d 4 xtr{[(U j dU ) 3 + (U dU j ) 3 } (AUA)} , (4.8) 

where, again, Stoke's theorem has been employed. As Uo is static we have Twz[^o] = and 
the remainder becomes a local object which is straightforwardly evaluated. 

4.1 Quantization of the three flavor collective Lagrangian 



In order to quantize the three flavor Lagrangian (|4.7|) we require the operators for spin and 
flavor as Noether charges. As a consequence of the hedgehog structure, the infinitesimal change 
under spatial rotations can be written as a derivative with respect to f2 

[rx9,[/(r,i)] = ^. (4.9) 

By the Noether construction this leads to the spin operator J = dL(A,Q a )/dfl. The quan- 
tization of the U SU(3) rigid top" proceeds by generalizing this result to the so-called right 
generators 



^ T -a 2 tl a = -J a , a=l,2,3 

R - = ( L Skyrme + L wz ) = { ~^ a , a=4,..,7 . (4.10) 

8Vt a N C B 



a=8 



The quantization prescription then demands the commutation relation [R a , Rb] = —ifabcRc 
with f a bc being the antisymmetric structure constants of SU(3). Explicit expressions for these 



generators in terms of an "Euler-angle" parameterization of A are presented in ref ||0[ . The so- 
called left generators, which are defined by the rotation L a = D^Rb, satisfy the commutation 
relations [L a , L b ] = if a b c L c . They provide the isospin, Jj = Li (i — 1, 2, 3) and hypercharge, 
Y = 2L 8 / v^3 operators. 

The generator R 8 is linearly connected to the so-called right hypercharge Y R = 2_R 8 /y / 3 = 1 
for B = 1 and Nc = 3. In analogy to the Gell-Mann Nishijima relation a right charge 

Qr = -Js + y (4.11) 

may be defined. Completing the analogy we note that the eigenvalues of Qr are 0, ±1/3, ±2/3, 
±1, .... Hence for Yr = 1 the relation (|4.11|) can only be fulfilled when the eigenvalue of J3 is 
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half-integer. This yields the important conclusion that the 577(3) model describes fermions. A 
priori this is not expected since the starting point has been an effective model of bosons. This 
discussion can be generalized to arbitrary Nc showing that the Skyrmion describes fermions 
when Nc is odd and bosons when N c is even [0O|1 . This, of course, is expected from considering 
baryons as being composed of N c quarks. We conclude that the proper incorporation of the 
anomaly structure of QCD leads to the desired spin-statistics relation. 

4.2 Flavor symmetry breaking and baryon spectrum 

For a realistic treatment of baryon states in the space of the collective coordinates we have 
to supplement the collective Lagrangian by the flavor symmetry breaking pieces associated with 
(|3.13| ). Substituting the flavor rotating hedgehog yields the symmetry breaking piece in the 
collective Lagrangian, 

L SB = (1 - D 88 ) (4.12) 

with the coefficient, 7 being linear in the symmetry breaking parameter 1 — x, i.e. 
32tt 



7 



-1) J dr {5'r 2 {l - cosF) - /3'cosF (f'V + 2sin 2 F) } . (4.13) 



D 88 (A) is defined in eq (|4.6|). Putting pieces (|4.7|) and ( [4.12|) together, the Hamiltonian for the 



collective coordinates is obtained as the Legendre transform H = — Yla=i Rafia — L 



H(A,R a ) = E d + l - 



1 1 

^2 ~ ]p 



1 ;C 2 (SU(3)) - -L + i 7 (1 - D m ) (4.14) 



2f3 2 y y " 8f3 2 2 



for B = 1 and Nc = 3. The constraint i?g = 2 ' wmcn yielded the spin-statistics relation, com- 
mutes with H permitting one to substitute this value. The term involving Y/ a =4 Ra h as been 
re-expressed by introducing the quadratic Casimir operator of SU(3), C 2 (5f/(3)) = X^ =1 72 2 . 
The standard SU{3) representations are eigenstates of C 2 (5?7(3)) with eigenvalues /i. For 
example, the octet representation 8 has fig = 3 while /iio = fijo = 6 and /z 2 7 = 8. These repre- 
sentations diagonalizef] the collective Hamiltonian in the absence of symmetry breaking, 7 = 0. 

Now consider the full collective Hamiltonian including the symmetry breaking. It seems 
reasonable to assume these 7 = eigenstates as a basis to diagonalize the full Hamiltonian. 
In a perturbative treatment up to third order in 7 for the | + baryons only the representations 



8, 10 and 27 contribute |42| . For that reason the perturbative treatment is still simple, although 
one must go beyond leading order. In particular this implies that the nucleon is no longer a 
pure octet state but rather contains sizable admixture of the nucleon type states in higher 
dimensional representations, 

\N) = |A^,8) + 0.07457 ( 5 2 |7V,TO)+0.04907 ( 5 2 |7V,27) + ... , (4.15) 

where the coefficients of the effective symmetry breaker 7/3 2 are computed from SU(3) Clebsch- 



Gordon coefficients | 43| . The nucleon is seen to have a roughly 25% amplitude to contain the 
10 state. 



1 The hedgehog structure of the classical configuration Uo constrains the permissible SU(3) irreducible rep- 
resentations to those which have at least one state with I = J [Ell . 
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Although this perturbative treatment provides a physical picture of the symmetry breaking 
effects it actually turns out that the full Hamiltonian ( |4.14j ) can be exactly diagonalized nu- 
merically. The important ingredient is that within a suitable "Euler-angle" representation of 
the rotations A, the symmetry breaker 1 — D 88 depends only on one of these eight angles. In 
each isospin channel the eigenvalue equation 

C 2 (SU(3)) + /3 2 7 (1 - D 88 )] * = esB ^ (4.16) 

then reduces to a set of coupled ordinary differential equations which can be integrated numer- 
ically. Here we do not wish to discuss this approach in full detail; rather we refer the reader to 
the original work by Yabu and Ando and exhaustive applications of this method involving 



the present authors |4Ti| . Having obtained the eigenvalue esB the baryon masses are 

straightforwardly computed from 



M B = E + - — -— j(j + l)-_ + — esB . (4.17) 





As already mentioned this diagonalization procedure is equivalent to the perturbation expan- 
sion. For small enough symmetry breaking /3 2 7 even first order is sufficient. In that (unjustified) 
case the famous Gell-Mann-Okubo mass formulae E7], RJ holds exactly: 



2(M N + M S ) = M E + 3M A (4.18) 
Mu - \h- = M H * - M s * = .1/v. - M A . (4.19) 

Additional corrections |55[ arise when we allow for non-zero classical .ff -meson fields to get 
induced by "rotations" Q a into the strange directions. These are energetically favorable since 
they maximize the strange moment of inertia (3 2 . With a parameterization 

the radial function W(r) is determined from applying a variational principle to j3 2 . In principle 
one must enforce that the ansatz ( [4.20| ) has no overlap with any global rotation of the classical 
solution (O). 



We adjust the only free parameter, e ~ 4 to the mass differences of the low-lying | + and 
| + baryons. The resulting baryon spectrum is shown in table |4.1| . Apparently the three flavor 
Skyrme model reasonably accounts for the empirical mass differences. The original studies 



48], |49], [50], [44j yielded far too low mass splittings between baryons of different strangeness 
for physically motivated parameters of the effective LagrangianQ. A major reason for the im- 
provement is the fact that 7 is significantly enlarged by including the effects associated with 
fx 7^ f-K [0- It is also apparent from table [4. X| that enforcing the zero overlap condition for the 
induced kaon components can be compensated by a small variation of the Skyrme parameter, 
e. This indicates that possible double counting effects play only a minor role. It is interesting 
to remark that the mass differences for the | + baryons deviate strongly from the predictions 
in leading order of the flavor symmetry breaking. This can easily be observed from the ratios 

(M A - M N ) : (M s - M A ) : (M 3 - M s ) = 1 : 0.52 : 0.85 (4.21) 

2 Many of these authors considered as a free parameter fitted to the absolute values of the baryon masses. 
Without the /?' term this yielded f n as low as 25MeV @. 
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Table 4.1: The mass differences, which are obtained by exact diagonalization of the collective Hamil- 
tonian (4.14), of the | + and | + baryons in the pseudoscalar model for e=4.o are compared to the 
experimental data. The values in parentheses are obtained by enforcing the zero overlap condition 
mentioned after (|4.2C| ) pOR . In that case the Skyrme parameter has slightly been readjusted to e=3.9. 
All data are in MeV. 



Baryons 


Model 


Expt. 


A 


- N 


154 (163) 


177 


£ 


- N 


242 (264) 


254 




- N 


366 (388) 


379 


A 


- N 


278 (268) 


293 


E* 


-N 


410 (406) 


446 




-N 


544 (545) 


591 


n 


- N 


677 (680) 


733 



which are in much better agreement with the experimental data (1:0.43:0.69) than the leading 
order result (1:1:0.5). Obviously the higher order contributions are important. This also indi- 
cates that the baryon wave-functions contain sizable admixture of higher dimensional SU (3) 
representations, cf. eq ( |4.15| ). Nevertheless the deviation from the Gell-Mann-Okubo relations 
( |4.18| ) is only moderate, in particular the equal spacing among the | + baryons is well repro- 
duced. Finally we note that, as discussed in point (iii) of section 2, the absolute mass of the 
nucleon is also too high in the three flavor case. Again we must rely on the N® corrections 
mentioned. 

4.3 Electromagnetic properties of baryons 

The value for the Skyrme parameter e = 4.0 obtained from this best fit to the baryon mass 
differences is next employed to evaluate static properties of baryons within this model. In order 
to do so one first constructs the Noether currents associated with the symmetry transformation 
( |3.9| ). A convenient method is to extend these global symmetries to local ones by introducing 
external gauge fields (e.g. the gauge fields of the electroweak interactions) into the total action 
i.e. (|1 . 1|) , ( [2.4|) , ( [3.12j ) and (|3.13|) . The Noether currents are then read off as the expressions 
which couple linearly to these gauge fields. This procedure is especially appropriate for the 
Wess-Zumino term ( |3.12|) because this non-local term can only be made gauge invariant by a 
lengthy iterative procedure j3"2"| . The final form of the nonet (a = 0, . . . , 8) vector (V®) and 
axial-vector (A®) currents reads [f|IJ (for Nc = 3) 



16tt 2 



c [ivpa 



tr {(tQ a e±?Q a t)pvP P P*} 



-ip'ti {Q a ({UM + MU\ c^} t- {MU + U^M,^}) } 



where Q a 



, 3 ) 2 ' • • • ' 2 ■ 



denote the Hermitian nonet generators. The combination 



Q c 



diag 



2 
3' 



1 1 

"3'~3 



(4.22) 



(4.23) 
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is of special interest because it enters the computation of the electromagnetic properties. The 
associated form factors of the \ + baryons (B) are defined by 



(B(p>)\Vr-\B(p)}=u(p>[ 



u(p) 



2M B 

Frequently it is convenient to introduce "electric" and "magnetic" form factors 

-,2 



q 



guq 2 ) = W) - 7^72 m* 



4M| 



GfM 2 ) 



Ft 



lit 



lit 



(4.24) 



(4.25) 



Substituting the rotating hedgehog configuration into the defining equation of the currents 
( |4.22| ) yields for the spatial components of the vector current^ 



V a 



\/3 

Vi{r)e ijk XjD ak + —B(r)e ijk fljX k D a8 + V 2 (r)e ijk x j d da pD aa Q l3 
jD 88 D ak + Vi{r)eij k x jddapDgaD a p + . . . , 



where 



B(r) 



2tt 2 



-1 F M 2 F 



(4.26) 



(4.27) 



is the baryon number density Q2.2| ). The explicit form of the radial functions Vi(r), . . . , V^(r) 
is given in appendix B of ref ||46|| . According to the quantization prescription ( |4.10| ) the an- 
gular velocities fl a are replaced by their expressions in terms of the right generators R a of 
SU(3). Taking the Fourier transform of the resulting matrix elements allows one to identify 
the magnetic form factor in the Breit frame ^] 



GfM 2 ) 



-8nM B / r 2 dr--j 1 (r\q\)\v 1 (r)(D e3 } B -—B(r)(D e8 R h 
Jo \q\ { 2cr 



h 



(4.28) 



-^V 2 (r)(d 3al3 D ea R ) B + V 3 {r)(D 88 D e3 ) B + V A {r)(d 3af3 D ea D w ) B | . 

Here the flavor index e refers to the "electromagnetic" direction (|4 . 2 3|) . The magnetic moment 
corresponds to the magnetic form factor at zero momentum transfer fi B = Gf 4 (0). Similarly the 
electric form factor is given by Fourier transforming the time component of the electromagnetic 
current 

G% = 4vr J o r 2 dr Jo (r\q\)l^-B(r)(D e3 ) B + — y 7 ( r )( J D ei i? i ) B + — V 8 {r){D ea R a ) B \. (4.29) 

The two new radial functions V-j(r) and V 8 (r) are listed in appendix B of ref as well. 
Integrating Vj and V 8 yields the moments of inertia, a 2 and (3 2 , respectively. Hence the electric 
charges are properly normalized. It should be remarked that the baryon matrix elements in 
the space of the collective coordinates are computed using the exact eigenstates of ( |4.16| ) and 
adopting the Euler-angle representations for the SU(3) generators |4(|. The results for the 
magnetic moments and the radii 



r 2 



6 dGUq 2 



fi B dq 2 



and r\ = —6 



dG B E {q 



dq 2 



(4.30) 



22 



Table 4.2: The electromagnetic properties of the baryons compared to the experimental data. The 
predictions of the Skyrme model are taken from ref j5| . 



B 


/x B (n.m.) 
e = 4.0 Expt. 


r 2 M (fm 2 ) 
e = 4.0 Expt. 


r|(fm 2 ) 
e = 4.0 Expt. 


P 


2.03 


2.79 


0.43 


0.74 


0.59 


0.74 


n 


-1.58 


-1.91 


0.46 


0.77 


-0.22 


-0.12 


A 


-0.71 


-0.61 


0.36 




-0.08 




£+ 


1.99 


2.42 


0.45 




0.59 




S° 


0.60 




0.36 




-0.02 




£- 


-0.79 


-1.16 


0.58 




-0.63 




E° 


-1.55 


-1.25 


0.38 




-0.15 






-0.64 


-0.69 


0.43 




-0.49 




S° -► A 


-1.39 


-1.61 


0.48 









are shown in table fLQ. As in the two flavor model |J the isovector part of the magnetic 
moments is underestimated while the isoscalar part is reasonably well reproduced. Despite the 
fact that the flavor symmetry breaking is large for the baryon wave-functions, the predicted 



magnetic moments do not strongly deviate from the SU(3) relations |53 



2 

2>A = -(/■«£+ + MS") = -2// S ° = Hn = fJ>S° = —^}iE°A ■ ( 4 -31) 

A more elaborate treatment of the flavor symmetry breaking is necessary in order to accommo- 
date the experimentally observed details of breaking the [/-spin symmetry which e.g. causes 
the approximate identity ~ fi p [p4]| . The moderate differences between the various mag- 



netic radii r 2 M is a further hint that symmetry breaking effects are mitigated. The comparison 
with the available empirical data for the radii shows that the predictions turn out too small in 
magnitude (except for the neutron electric radius). This is a strong indication that essential 
ingredients are still missing in the model. In section 5 it will be explained that the effects, 
which are associated with vector meson dominance (VMD), will account for this deficiency. 
Nevertheless the overall picture gained for the electromagnetic properties of the ^ and | 
can at least be characterized as satisfactory, especially in view of the fact that the only free 
parameter of the model has been fixed beforehand. 

4.4 Effects of symmetry breaking on baryon matrix elements 
and strangeness in the nucleon 

Theorists typically wish for symmetry breaking effects to be negligible (the so-called "spher- 
ical cow" approximation) but Nature says otherwise. This is very apparent in the case of low 
energy strong interactions (QCD). The Gell-Mann-Okubo mass formulae, which amount to ap- 
plications of a Wigner-Eckart theorem for first order \$ type symmetry breaking, furnish sum 
rules rather than a complete description. As we discussed in section 4.2, the Skyrme model 
provides a non-trivial playground for treating symmetry breaking. The Yabu-Ando equation 

3 The conventions are i, j, k — 1, 2, 3 and a, (3 = 4, . . . , 7. 
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( 4.16Q gives an exact (within the model) wave-function for each baryon state at any strength 
of the symmetry breaking parameter 7/3 2 (oc underlying quark masses). The physical results 
vary smoothly with 7/3 2 (although higher quantum corrections would be expected to give weak 
non-analytic corrections). 

The physical interpretation of symmetry breaking in the model may be seen from eq fl4.15| ). 
The higher SU(3) representation components in the baryon wave-function can only emerge in 
a quark framework by having quark-ant iquark pairs present in addition to the three "valence" 
quarks. Clearly such effects would be difficult to treat in the non-relativistic quark model 
approach. On the other hand, it should be recognized that the Skyrme model approach is 
based on a collective semi-classical treatment. 

In the last few years there has been a greatly renewed interest in the study of symmetry 
breaking effects for ordinary nucleons. This was stimulated by new experiments on polarized 



lepton deep inelastic scattering off nucleons |55| which seem to indicate that the pure valence 
quark picture of the nucleon has serious drawbacks. The Skyrme model has the advantage of 
giving a simple and roughly accurate quantitative explanation of these experiments. In detail 
one needs the strangeness conserving proton matrix elements 



[P(p)\q l l,l 5 qi\P(p))=u(p') 



2 1 ir«^ 



lr,u (p) 



(4.32) 



for this discussion. Of related interest are the flavor changing matrix elements 

(B'(p')\V«\B(p)) = u(p')[ ltl g v (q 2 ) + ..]u(p) , 
(B'(pi)\A«\B(p)) = u(p>) \ W g A (q 2 ) + ..] u{p) . 



(4.33) 



between different baryons (B',B). Here we have omitted contributions proportional to the 
momentum transfer q^. 

Knowledge of the g A (B, B') and gv{B, B') are crucial for the theory of baryon semi-leptonic 
decays. First let us consider the calculation of the axial vector matrix elements g A (B, B'). Our 
main interest in this brief discussion will be to examine the effects of symmetry breaking. The 
leading order term (in 1/Nc) of the spatial components of the axial current is straightforwardly 
obtained to be 



J d 3 rA« = CD ai (A) 



(4.34) 



where A(t) is the collective coordinate matrix. The constant C denotes an integral over the 
chiral angle. We refer the interested reader to refs |57], for the explicit expression. Then, 



g A (B',B)=C(B'\D a3 \B) 



(4.35) 



The flavor index a has to be chosen according to whether strangeness conserving (a = 1, 2, 3, 8) 
or strangeness changing (a = A, ... ,7) processes are being considered. The corresponding 
result for the axial charge of the nucleon g A = g\ +l2 (p, n )i as measured in neutron beta-decay, 
is predicted too low in many soliton models. This problem is already encountered in the two 
flavor model and gets worse in SU (3) as the Clebsch-Gordon coefficient associated with Di +i2 3 
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Table 4.3: The matrix elements of the axial-vector current ( [4.35 ) between different baryon states 
in the flavor symmetric limit. Displayed are both the strangeness conserving (a) and strangeness 
changing (b) processes. The first column gives the relevant flavor component of the axial current. 



A 7 " 


n — > p 
F + D 


£- -> A 


(a) 

£- -E° 


D-F 




A — > p 
^F + D) 


£~ — > n 
D-F 


(b) - A 


E~ — > S° 
75^ + ^) 



changes by a factor of 7/10. As symmetry breaking is increased the SU(3) prediction for qa 
becomes larger [42 

g A {SU{3)) = ^ [l + 0.0514 7 /5 2 + . . .] g A (SU(2)) . (4.36) 

Actually the exact treatment shows that with increasing symmetry breaking the two flavor 
result is approached, although only slowly. Taking everything together, including subleading 
terms in ( |4.35|) , finally gives g A = 0.98 for e = 4.0 which is about 4/5 of the experimental 
value g A (expt.) = 1.26. 

We may understand the tendency to approach the SU(2) limit for large 7/? 2 as follows. In 
the small SU(3) breaking case, there is just a small extra "cost" for producing an ss pair rather 
then a uu or dd pair. As 7/3 2 gets larger it is more expensive to make an ss pair and eventually 
ss pairs should be absent from the nucleon wave-function, recovering the SU(2) picture. 

Returning to the general case one should first note that flavor symmetry relates the octet 
axial current matrix elements between various baryons. Conventionally they are expressed 
using SU(3) covariance in terms of two unknown constants (or reduced matrix elements) F 
and D. One has to use models to determine these constants. In the flavor symmetric Skyrme 
model one finds |BB[ D/F = 9/5 and D + F = 7C/15 = g^. In table [4.3| the flavor symmetric 



dependences of the axial matrix elements on F and D are displayed. As one departs from the 
flavor symmetric case the baryon wave-functions acquire admixture from higher dimensional 
SU(3) representations making the SU(3) covariant parameterization in terms of F and D 
inadequate. In the presence of SU(3) symmetry breaking we must, without further assumptions, 
parameterize each decay amplitude separately. It is still reasonable to maintain the isospin 
invariance relations. 

As an example of the perturbative corrections consider the axial A — * p transition in the 
Cabibbo scheme for semi-leptonic hyperon decays. The analog of ( |4.15|) for the A hyperon 



is 

|A) = |A,8) + |-7/3 2 |A,27) + ... . (4.37) 
50 

Noting that the D-functions mix different SU (3) representations, we get 

(p T \D K - 3 \A T) = ^= - ^| 7 /? 2 + • • ■ , D K - 3 = -L (^3 - zD 53 ) . (4.38) 
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B' I A°, I B 




Figure 4. 1: The variation of axial 
vector matrix elements with the 
effective symmetry breaking pa- 
rameter 7/3 2 . 

Full line: (p|s7 3 7 5 s|p); dashed 
dotted line: (p|w7375s|A); dotted 
line: (n|w7375s|£~); long dashed 
line: (A|w7 3 7 5 s|S~); dashed line: 

{jp\ujzlhd\ n ) ■ These matrix ele- 
ments, which are taken from refs 
DjJ and 



581, are normalized to 



the flavor symmetric values. 



7-/S 2 



Of course, this expansion just provides a first approximation to the symmetry breaking depen- 
dence of the Cabibbo matrix elements. Using the exact treatment initiated by Yabu and Ando 



44] this dependence can be computed numerically as shown in figure 14 for some processes 
of interest |)7|, BS|. Those results are normalized to the SU(3) symmetric values in table fO 
to illustrate that the matrix elements vary in different ways with symmetry breaking. In this 
figure also the variation of the nucleon matrix element of the flavor conserving axial current 
^3(0) = (A^|s7 3 7 5 s|A^) is displayed. Obviously H 3 (0) decreases very rapidly with increasing 
symmetry breaking. This is easily visualized, as mentioned above, as a reflection of the in- 
creased cost of making extra ss pairs in the nucleon wave-function as 7/3 2 increases. On the 
contrary the Cabibbo matrix elements exhibit only a moderate dependence on 7/3 2 . It is this 
different behavior of the matrix elements that makes the application of exact flavor symme- 
try to the analyses of the EMC-SLAC-SMC experiments suspicious. Stated otherwise, the 
strange quark contribution to the nucleon matrix element of the axial singlet current (loosely 
"proton spin" ) may be decreased significantly as a consequence of symmetry breaking without 
contradicting the successful Cabibbo scheme for the semi-leptonic decays of the hyperons. 

As an interesting contrast to the axial matrix elements, consider the evaluation of the 
vector matrix elements gv{B, B') needed for the hyperon semi-leptonic decays. The dominant 
contribution is given by the matrix elements of the SU(3) flavor generators 

g«(B',B) = (B'\L a \B) . (4.39) 

For example if we sandwich the generators L^- between the perturbative A and p states given in 
(|4.37| ) and (|4.45| ) and recognize that group generators can only connect states belonging to the 
same irreducible representation, we see that symmetry breaking corrections start out as (7/9 2 ) 2 
rather than 7/? 2 . This is just a demonstration of the Ademollo-Gatto theorem [S3], which 



"protects" the vector matrix elements against small symmetry breaking corrections. Since 7/3 2 
is large the numerical validity of this result is questionable. However, the exact Yabu-Ando 
scheme does confirm |57j that vector matrix elements suffer at most, 40% deviation from the 
symmetric values, even for large symmetry breaking, e.g. 7/5 2 ~ 7. 

A reduction of strangeness in the nucleon with increasing 7/3 2 is also predicted for the scalar 
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strange content fraction of the proton 



X. 



(p\ss\p) - (0\ss\0) 



(p\uu + dd + ss\p) — (0\uu + dd + ss\0) 



(4.40) 



Here the state |0) refers to the soliton being absent. Models of quark flavor dynamics, as e.g 
the one of Nambu-Jona-Lasinio [[J, indicate that matrix elements of quark bilinears q\ a q may 
be taken as proportional to the matrix elements of tr 
wardly get 



A a [U + — 2) . Then we straightfor- 



X x 



43 
2250 



7/? 2 + ... • 



(4.41) 



In this case, however, the deviation from the flavor symmetric result |K| (X s = 7/30) is 
considerably mitigated as compared to the variation of H s defined in ( [4.32|) . The symmetry 
breaking has to be as large as 7/? 2 ~ 4.5 to obtain a reduction of the order of 50%. In the case 
of H s this was already achieved for 7/3 2 ^2.5. In any event, the additional quark-ant iquark 
excitations in the nucleon, which are parametrized by the admixture of higher dimensional 
SU(3) representations (^.15|) , clearly tend to cancel the virtual strange quarks of the pure octet 
nucleon. 

The three flavor Skyrme model under present consideration provides a convenient way to 
study the nucleon matrix elements of the vector current sj^s . These are theoretically inter- 
esting because they would vanish in a pure valence quark model of the nucleon and so test 
finer details of nucleon structure. They are experimentally interesting because they can be 
extracted from measurements of the parity violating asymmetry in the elastic scattering of 
polarized electrons from the proton. The precise form factors needed are defined by 



(P(p')\q s l,qs\P(p))=u(p' 



l»F s (q 2 ) + 



2M r 



-F s (q 2 



u{P) 



(4.42) 



These form factors are currently under intensive experimental investigation, cf. refs and 
have been estimated in various models. The models range from vector-meson-pole fits |)3| 
of dispersion relations [53 through vector meson dominance approaches and kaon-loop 



calculations with [65] and without [^6] vector meson dominance contributions to soliton model 
calculations |B7|, BSfl. The numerical results for the strange magnetic moment fi$ = F s (0) w 
—0.31 ± 0.09 . . . 0.25 are quite diverse. The predictions for the strange charge radius r| = 
— 6dF s (q 2 ) / dq 2 \ q=0 are almost equally scattered r| ~ —0.20 . . . 0.14fm 2 . 

In order to evaluate these form factors in the three flavor Skyrme model one requires the 
matrix elements of the "strange" combination 



Q' 



(4.43) 



rather than the electromagnetic one ( [4.23| ) between proton states. Using the same value e = 4.0 
as used consistently for the three flavor pseudoscalar model yields the predictions 



Us = — 0.13n.m. 



-O.lOfm 2 



(4.44) 



Here n.m. stands for nuclear magnetons. It should be stressed that these results are obtained 
within the Yabu-Ando approach, i.e. the proton wave-function contains sizable admixture of 
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higher dimensional representations. If a pure octet wave-function were employed to compute 
the matrix elements of the collective operators the strange magnetic moment would have been 
Us = —0.33. The proper inclusion of symmetry breaking into the nucleon wave-function is 
again seen to reduce the effect of the strange degrees of freedom in the nucleon. We already 
discussed above that as the strange quarks within the nucleon become more massive (effect of 
symmetry breaking) their excitation becomes less likely. 

In the next few years a great deal of new experimental information on the form factors F s 
and F s should become available. This would enable more accurate comparison with (for a given 
effective meson Lagrangian) the essentially parameter free predictions of the soliton theory. 

5. The nucleon as a vector meson soliton 

According to the modern view of the Skyrme model approach we should start from the 
"full" effective Lagrangian which contains mesons of all spins. The practical criterion on which 
particles to include is to find an effective Lagrangian which does a good job of explaining the 
low energy experimental data in the meson sector. On these grounds it is evident that the 
vector mesons should be included in the effective Lagrangian. In this section we give a brief 
sketch of the soliton sector of the Lagrangian of pseudoscalars and vectors (see section 3.3) and 
note that it leads to significant improvements of many predictions. In particular it is crucial 
for discussing the so-called proton spin puzzle. 

5.1 Generalized soliton ansatz and profile functions 

As a first step we construct the soliton of the Lagrangian defined in section 3.3. The gener- 
alization of the hedgehog ansatz (|2.6|) to the vector meson model requires the time component 
of the uj field and the space components of the p field to be different from zero. Parity and 
grand spin symmetry^] allow for three radial functions 

& = exp (If ■ rF(r)) , = ^ , p\ = ^-e ija rj . (5.1) 

Substituting these ansatze into the action described in section 3.3 yields the classical mass, 

-f 2 r 2 1 G 2 



E 



4vr / dr[f(F' 2 r 2 + 2sin 2 F) - ^{J 2 + m 2 ^ 2 ) + \[G' 2 + ^(G + 2) 2 ] 



2g 2 ' p ' g 2 

+ ^f(l +G- cosF) 2 + ^F'usm 2 F - ^G'usmF 
9 2 9 9 

+—F'uG(G + 2) + -( 72 + -fy)F'uo[l - 2(G + l)cosF + cos 2 F] 

r a' 
+ (1 - cosF){45V 2 + 2(2/3' - -^)(F' 2 r 2 + 2sin 2 F) 

9 

Icy' > i 

-\co 2 r 2 -2(G+1- cosF) 2 - 4(1 + cosF)(l + G — cosF)]} . (5.2) 

9 

Application of the variational principle to this functional yields second order coupled non-linear 
differential equations for the radial functions F(r), u(r) and G(r). The boundary conditions 

1 The grand spin J +1 is the characteristic invariance of the hedgehog ansatz. 



28 



for the chiral angle F(r) = it and F(oo) = 0, which correspond to unit baryon number, also 
determine the boundary conditions of the vector meson profiles via the differential equations 
and the requirement of finite energy. For example we find G(0) = —2. A typical set of resulting 
profile functions is shown in figure [5.1| . 




Figure 5.1: The profile func- 
tions which minimize the clas- 
sical vector meson energy func- 
tional (|5.2j) for the parameters 
g = 5.85, h = 0.4, g vv $ = 1.9 
and k = 1.0. u(r) is measured 
in units of m p . 



As in the pseudoscalar model we have to generate states with good spin and isospin from 
this classical field configuration. To start with, one introduces collective coordinates Q2.8Q for 
all fields which have non-vanishing spin or isospin. However, an additional complication arises 
because there are vector meson field components which vanish classically but get excited by 
the collective rotation. In the two flavor case the appropriate ansatz for these excitations reads 



1 



Po = -A{t) [fc(r)n + 6(r)(f ■ n)f] A\t) 



<3>(r) 

~2g~ 



(5.3) 



where the angular velocity of the rotating soliton, fij is defined in ( |2.10| ). The three radial 
functions 6>6 and $ are not the only ones which get excited. As these radial functions are 
non-zero they provide sources for the non-strange component of the iso-singlet pseudoscalar 
field via the e-terms ( 3.29Q . In the two flavor formulation the appropriate ansatz which takes 
into account the pseudoscalar nature of the rj field reads 



U{r,t) 



Jvt (r) 



A(t)U (r)A j (t) with r) T (r) = —n(r)r ■ ft 

Jit 



(5.4) 



As we will observe shortly, the excitation of this r] field plays a decisive role in the context of 
the proton spin puzzle. The additional radial functions are determined from extremizing the 
moment of inertia for rotations in coordinate space, 



(fi + 6) 2 + 2(6 - 1 + cosF) 2 



(3^1 + 2££ + & )r 2 + 2(G 2 + 2G + 2)& + AG\^ + 66 - 26 - £ 2 + 1 




+ -7i0F'sin 2 F + -^F' (G - 6)(1 - cosF) + (1 - cosF) 2 - 
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'smF(G - 6 + 2 - 2cosF) + 0sinF(^ - G') 



+(j)F l [2 + 2sin 2 F + (6 - G - 2)cosF - 2(6 + 6)] } 
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L {V [(G + &)G + + - cosF) 2 - 2GcosF]] + 77 (G^ - G'^)} j 



(5.5) 

together with suitable boundary conditions. In eq ( |5.5| ) we have not displayed the explicit 
contributions from the symmetry breakers (which are in fact small). We will mostly limit the 
present discussion to the two flavor case. In the case of three flavor vector meson models the sit- 
uation is even more complicated as also K* type fields get excited. Also there will be additional 
symmetry breakers on the level of the collective Lagrangian which are of the form £< =1 L> 8 A 
and stem from terms which are linear in the time derivative. They can straightforwardly be 
implemented in the collective quantization approach. Here we will omit details but rather refer 
the reader to the literature p7| , |4T3| . The general pattern for computing baryon properties is 
essentially the same as that discussed for the Lagrangian of only pseudoscalars in section 4. 

5.2 Axial singlet current and proton spin puzzle 

Notice that in ( |5.5| ) we included by hand a mass term for the rotationally excited profile 
77 (r) of a pseudoscalar isosinglet field. Actually the existence of such a term has not yet been 
justified. Before proceeding we must do so since the term turns out to be very important. 

In section 3.1 we mentioned that the QCD axial singlet current 

J ti = m7m75« + dl^d + s7m75S, (5.6) 

is not conserved even for zero quark masses: d^J^ = G, where the £74(1) anomaly G is 
proportional to the product of the QCD field strength tensor and its dual. In order to 



mock up this non-conservation equation at the effective Lagrangian level |69[ we may add the 
terms 

c^ 2 *G n / det U \ 

2 G + T2 ta (suJtj' (57) 

where G is now considered a composite glueball field which "dominates" the {7a (1) anomaly. 
Here we assumed three light flavors and also that the strong CP violation parameter 9 is zero. 
Furthermore it is necessary that, except for the terms representing quark mass symmetry break- 
ing, all the other terms in the effective Lagrangian be invariant under {7a(1). The parameter c 
above is determined by 

m i « e^*. (5 - 8) 

in the approximation where the quark mass terms are neglected. 770 is the SU(3) singlet 
pseudoscalar field as in (|4.1|) . This equation arises after noting that G is like an auxiliary field 
and may be integrated out: G = i] / (^/Ecf^) . In the effective Lagrangian the realization of 
the axial singlet current, obtained by a Noether variation, is 

4°l = Veud.vo + j® . (5.9) 
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Here the first term is the contribution from the pseudoscalar field and the second term is due 

7(0) 

'5,n 



to the addition of vector fields. has a complicated structure but, in particular, contains 



non-derivative terms like e^ ual3 ti{p v p a pp). Using this decomposition we may write the equation 
of motion for the r/ field as 

(d 2 + m^rjo = -^-^Jg, (5.10) 

which shows that the vector meson contribution to the axial singlet current may act as a source 
for a non-trivial excitation associated with the rjo field in the soliton sector. 

Now the form factors for the proton matrix elements of the axial singlet current are obviously 
just the sums of the three separate form factors introduced in ( |4.32|) : 



H(q 2 )=H u (q 2 ) + H d (q 2 ) + H s (q 2 ) and H{q 2 ) = H u {q 2 ) + H d (q 2 ) + H s {q 2 ). (5.11) 



If the vector mesons are not present, eq (|5.9|) shows that the operator for the axial singlet 
current must be (even in the soliton sector) a pure derivative. This means that, regardless 
of the details of the calculation, the matrix element for the sum of the three terms in ( |4.32| ) 



must be proportional to the momentum transfer q^ . Thus H(q 2 ) is non-zero and H(q 2 ) = 0. 
From the theory of Dirac particles we recognize that the quantity H(0) has the interpretation of 
twice the quark spin part of the proton's angular momentum. We see that the Skyrme model of 
only pseudoscalars predicts that the expectation value of the net quark spin operator vanishes; 
the total angular momentum (1/2) of the proton must involve, at a fundamental level, the 
rotational and gluonic pieces! Note that the above argument for H(0) = with the Lagrangian 
of only pseudoscalars continues to hold even if symmetry breaking contributions are taken into 
account f58fl . 

The situation is a little different when vector mesons are included in the effective Lagrangian. 
Since has pieces which are not pure derivatives it then is possible to obtain H(0) 7^ . A 
convenient parameterization for this calculation in the effective Lagrangian model is 

W)\jf}\P{p)) = Hiq 2 )^) . (5.12) 

Once all the radial functions have been determined as before from the appropriate variational 
principles, it is straightforward to compute H(0) from eq ( |5.12| ). One only has to recall that 



under the collective coordinate quantization the angular momentum operator is given by J = 
fl/a 2 . The numerical results for a variety of allowed parameters (cf. discussion after eq ( |3.29| )) 
are displayed in table f).l\ Surprisingly the predictions of the vector meson model for H(0) are 

Table 5.1: Predictions for the matrix element of the proton axial singlet current for various allowed 
sets of parameters in the vector meson model. Results are taken from ref 



h 


0.4 


0.4 


0.4 


0.7 


0.5 


0.2 


0.1 


gvv® 


1.9 


1.9 


1.9 


2.2 


2.0 


1.7 


1.5 


K 


0.0 


0.5 


1.0 


0.0 


0.0 


0.0 


0.0 


F(0) 


0.34 


0.33 


0.30 


0.29 


0.34 


0.32 


0.28 



very robust against possible changes of the parameters of the model. 
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Even though we get a non-zero value for H(0) in the vector meson model it is still small 
compared to H(0) = 1, the expectation from the simple non-relativistic quark model. Quali- 
tatively the soliton model results with and without vectors are similar. Since one has a natural 
prejudice that the quark model results should be roughly correct, this would at first seem to 
be a serious defect of the soliton approach to nucleon properties. 

Of course one can only make an accurate judgment on the matter by appealing to experi- 
ment. H(0) can be found from eq ( |5.11| ) if we can experimentally obtain separately H u (0), H d (0) 



and H s (0). The linear combination 

H u {0) - H d {0) = g A = 1.257, (5.13) 

is reliably obtained by an isotopic spin rotation of the axial form factor describing neutron 
beta-decay. Similarly the estimate for the "eighth" octet component 

H u (0) + H d (0) - 2F S (0) « 0.575 ± 0.016, (5.14) 

may be gotten from an SU(3) flavor rotation of the data on hyperon beta-decay experiments. 
Clearly one more linear combination is needed in order to disentangle the individual i?j(0) 
and that situation existed for many years. About ten years ago different experimental groups 
(EMC, SLAC, SMC) used polarized lepton beams to probe the structure of nucleons. The deep 
inelastic scattering data |55j were used to extract the combination 



4H U (0)+H d (0)+H a (0), (5.15) 

in which the axial current form factor for each quark is weighted proportionally to the square of 
the quark electric charge. Combining these data resulted in H(0) ~ 0.3 [[71]. The experimental 



results were later on confirmed 1721, 1731, 1741. However, it turns out that the theoretical extraction 



of H(0) is quite complicated as it involves a careful treatment of perturbative QCD corrections. 
The value J75J 



H(0) = 0.27 ±0.04 (5.16) 

is nowadays considered correct. At the time this low value was considered hard to understand 
and the situation was called the proton spin puzzle. We have just seen that the soliton approach 
does however provide a simple explanation of such a low value. 

Clearly the prediction of the vector meson treatment described in Table 5.1, yielding H(0) 
about 0.30, is in good agreement with the data. From this we learn two things. First, the sim- 
plest quark model does not give a good description of the spin structure of the nucleon. Second, 
the soliton approach based on an effective Lagrangian including vector mesons markedly im- 
proves the qualitatively reasonable predictions of the soliton treatment based on a pseudoscalars 
only effective Lagrangian. A physical interpretation of the latter statement is that the pseu- 
doscalars only Lagrangian mainly probes the "pion cloud" of the nucleon while the vector 
Lagrangian probes a little more deeply. 

For completeness we remark on a possible caveat. The estimate of (|5.14j) is based on the 
use of exact SU(3) symmetry. However in Fig 4.1 of section 4.4 we showed that precisely this 
current matrix element is expected to exhibit stronger suppression than others due to SU(3) 
symmetry breaking. Nevertheless it turns out that [170] the numerical evaluation of -f^(0) is 
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not very sensitive to this feature. This is to be contrasted with the behavior of H s (0), which 
decreases rapidly with symmetry breaking, cf. section 4.4. 



5.3 Other improvements with vector mesons 

The famous problem of explaining the neutron-proton mass difference is another one which 
requires the addition of vector mesons to the effective Lagrangian in order to obtain a sat- 
isfactory solution in the nucleon-as-soliton picture. It is known that the electromagnetic 
(i.e. one photon loop) contribution has the wrong sign. After correcting for the electromag- 
netic interaction the remaining "strong" part of the neutron-proton mass-difference should be 
(M n — M p ) strong ps (2.0 ± 0.3)MeV [f76| . At the quark level this arises from the down quark-up 
quark mass difference, controlled by the parameter y in eqs ( |3.6|) and ( |3.7|) . Information on y 
can be most easily gained by analyzing the K + -K° mass-difference, yielding y ps (—0.4... — 0.2) 
33| . To understand the problem it is helpful to consider the contribution of the (presumably 



dominating) 5'-type symmetry breaker to the neutron-proton mass-difference. Since the d—u 
quark mass difference clearly exists with only two flavors it is interesting to first consider the 
problem at this level. Then the relevant piece of the 5' term is proportional to 

tr [r 3 (U + C/ f )] . (5.17) 

Using the ansatz (|5.4j ) we see that U = exp(it]T)[cos(ilj) + in ■ rsin(^)], where ip is some angle. 
Then ( |5.17| ) is proportional to sin(?7r). In other words the contribution vanishes unless the 
field t]t gets excited due to the collective rotation (or any kind of symmetry breaking). Now 
( 5.1 U| ) together with ( |5.9j ) shows that this will not happen if only pseudoscalars are present 



in the effective Lagrangian; the vector meson contribution must also be present. This is 
analogous to the situation concerning the proton spin puzzle. The contribution of the 6' term 
turns out to be 

(M n - M p ) strong = J d 3 r smF(r) V (r) + ... . (5.18) 

Using the full two-flavor vector meson result for which was already employed to compute 
H(0) yields 



(M n - M p ) strong ps 1.4 MeV (5.19) 

which, not surprisingly, turns out to be about as robust against changes of the parameters as 
is H(0). This prediction is still somewhat too small when compared to the empirical value. 
However, it turns out that the missing ~ 0.5MeV can be attributed to three flavor effects as 
matrix elements of D 38 are non-vanishingQ 

The addition of vector mesons also plays an important role in the discussion of the "sizes" 



of the nucleons: the nucleon radii. As can be observed from table [4.2| the Skyrme model of 



pseudoscalars only seriously underestimates the empirical values for the baryon radii. The 



presence of the to meson provides an increase of the isoscalar radius |52] 



(r 2 ) /=0 ~ (r 2 ) B + A , (5.20) 
mi 



2 



It should be remarked that {p\D^\p) quickly approaches zero as SU(3) symmetry breaking is increased. 



This decreases ST/ (3) type contributions to (5. IE) 
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where (r )b is the radius associated with baryon number current ( [2.2|) . The additional piece in 



eq ( |5.20|) is a consequence of (approximate) vector meson dominance in this model [78|, which 



indeed is observed when including the vector mesons in a chirally invariant manner. As can be 
seen from table [4.2| this increase of about 0.35fm 2 will significantly improve the predictions for 



the radii. 

A similar interesting improvement due to vector mesons is obtained in the context of meson- 
baryon scattering. In these investigations one introduces small fluctuations off the classical 
soliton. Eventually these fluctuations are quantized to represent in- and out-going meson 
fields, thereby determining the scattering matrix J75| . It turns out that in the pseudoscalar 
Skyrme model the phase-shifts extracted from this scattering matrix rise almost linearly with 
the momentum of the in-going pion. This undesired feature is mostly due to the contact 
interaction between pions contained in the Skyrme model Lagrangian (cf. section 2). When 
introducing vector mesons this contact interaction is essentially replaced by the exchange of 
such a vector meson, 

— 2-^^-2- ( 5 - 21 ) 

m z p q z — m z p 

As this interaction decreases for large momentum transfers, q 2 the resulting phase-shifts assume 
a constant value for large energies rather than rising linearly [BD[ . Clearly this effect is similar to 



the one observed when going from the Fermi to the standard model of electro-weak interactions. 

These examples show that while the inclusion of vector meson degrees of freedom involves 
quite a few technical details it clearly provides a more realistic picture of the nucleon as a chiral 
soliton. 

6. Summary and discussion 

Aside from the mass spectra and current matrix elements of the low-lying | + and | + baryons 
treated here the soliton approach has been extensively employed to study meson nucleon scat- 
tering [79 , jSOfl , baryons containing a heavy quark nucleon-nucleon scattering [3^|, few 



nucleon systems |83| and nuclear matter 

In the present survey, we started out with a historical introduction (section 1) and a concise 
technical summary of the original two flavor Skyrme model (section 2). In these sections the 
physical interpretation and justification of the model were emphasized; it is hoped that they 
will be useful to beginners in this area of research (see also and 

We next attempted to develop the generalization of the original Skyrme model which is 
suggested by the large Nq approximation to QCD. In this approach the Skyrme Lagrangian is 
to be replaced by a more general effective Lagrangian containing mesons of all spins. Perhaps 
some day an analytic expression in this framework will be found. At present it seems necessary 
to obtain an approximation based on including the lowest energy resonances and constraining 
the model by the symmetries of the underlying QCD. The concept of chiral symmetry which 
plays a crucial role in this extension was explored in section 3. Furthermore the original Skyrme 
model of two light flavors was extended to three flavors, as it is now well established that the 
nucleons belong to a flavor SU(3) multiplet. 

It is worthwhile to stress that once the effective Lagrangian has been determined from the 
meson sector, the soliton approach provides in principle a zero parameter description of baryon 
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properties (In our case we introduced just one parameter which had to be fit from the baryon 
sector.). 

In section 4 we studied the technical tools needed to treat the flavor SU(3) symmetry 
and its breaking at the (collective) baryon level. These were applied to the calculation of 
various interesting baryon matrix elements. Finally section 5 sketched the treatment of baryons 
based on an effective Lagrangian which also included the vector mesons. An application to 
the so-called proton spin puzzle demonstrated that the soliton approach seems to give a neat 
description of, otherwise hard to explain, experimental results on the quark spin structure of 
the nucleon. The improvements one encounters on including the vector mesons are in accord 
with the intuitive notion that the addition of higher mass resonances in the meson sector leads 
to a progressively more detailed understanding of the short distance structure of the nucleon- 
as-soliton. 

We are happy to acknowledge the stimulating interactions we have had with many collabo- 
rators and colleagues while doing research related to the topics reviewed here. 
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